function [params, paramsParts] = load_param_estimates(data, indivRunFolders)
	%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	% This function loads the parameter estimates from the 3 models + combine them into an object called params.
	% It also splits them by parts and combine that into an object called paramsParts.
	%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	%%%%% Inputs:
	% data:								object (necessary to get dims that are used to split params by parts)
	%	providerAdoptions:						object
	%		.dims
	%	providerExits:							object
	%		.dims
	%	consumerAdoptions:					object
	%		.dims
	%
	% indivRunFolders:					object:
	%	.consumerAdoptions:					string
	%	.providerAdoptions:					string
	%	.providerExits:						string
	%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	%%%%% Outputs:
	% params:							object:
	%	.consumerAdoptions:					NumParams x 1
	%	.providerAdoptions:					NumParams x 1
	%	.providerExits:						NumParams x 1
	% paramsParts:						object:
	%	.providerAdoptions:					cell(NumParts,1)
	%	.providerExits:						cell(NumParts,1)
	%	.consumerAdoptions:					cell(NumParts,1)
	%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	
	%%%%%%%%%%%%%%% Load parameter estimates from the 3 models + split by parts and save into object params  %%%%%%%%%%%%%%%
	%%%  Load parameter estimates from the 3 submodels + put them into object params
	load(sprintf('%s/results.mat', indivRunFolders.providerAdoptions), 'params_hat');
	params.providerAdoptions = params_hat;
	clear params_hat;
	
	load(sprintf('%s/results.mat', indivRunFolders.providerExits), 'params_hat');
	params.providerExits = params_hat;
	clear params_hat;
	
	load(sprintf('%s/results.mat', indivRunFolders.consumerAdoptions), 'params_hat');
	params.consumerAdoptions = params_hat;
	clear params_hat;
	
	%%% Split params_parts for the 3 submodels + save everything in object paramsParts
	paramsParts.providerAdoptions  = split_combine_parts(data.providerAdoptions.dims,  1, params.providerAdoptions);
	paramsParts.providerExits      = split_combine_parts(data.providerExits.dims,      1, params.providerExits);
	paramsParts.consumerAdoptions = split_combine_parts(data.consumerAdoptions.dims, 1, params.consumerAdoptions);
	
	clear providerAdoptions_params providerExits_params consumerAdoptions_params;
end
